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Abstract 

The dielectric response in relaxor ferroelectrics is analyzed in the framework 
a model for the polarization dynamics in the presence of polar clusters. We as- 



■ sociate the origin of polar clusters with the atoms displaced from their centro- 

■ symmerical positions even above Tc. Their collective hopping in multi-well 
potentials induced by disorder is analogous with the situation in glasses. The 

l> . 

0^ . theory explicitly takes into account the distribution of cluster reorientation 



frequencies and the effect of cluster-cluster interactions in highly polarizable 



. crystals, which we describe in terms of the local field distribution function. 

O ■ The dielectric constant is obtained from an integral master equation for the 



polarization dynamics in the presence of a time dependent electric field. The 
theory is applied for the analysis of the shape of the frequency dependent 
permittivity in the typical relaxor ferrolectrics PST as a function of tempera- 
ture. The comparison of the theory with experiment shows that in contrast to 
earlier assumptions, the observed Vogel-Pulcher dependence of the permittiv- 
ity maximum is a consequence of the Vogel-Fulcher temperature dependence 
of the cluster reorientation frequency. 

PACS numbers: 77.80.-e, 64.60.-i 



1 



Typeset using REVT^ 



I. INTRODUCTION 



In the family of ferroelectrics a large group of mixed and disordered "relaxors" has 
been identified including perovskites PhMgii^Nh2/2,0'i (PMN), PbSci/2Tai/203 (PST), La- 
modified PhZri-xTixOz (PLZT), tungsten-bronze structure oxides like PhxBai-xNh20Q 
(PBN), etc. P|^T3| Due to the effect of configurational disorder the properties of relaxor 
ferroelectrics are very different from those in translationally invariant ferroelectrics. Among 
the unusual properties of these materials are a) the coexistence of slow kinetics typical for 
spin glasses with a very large dielectric constant indicating intermideate range polar order 
on the nanometer length scale which can be transformed to true long range order, {i.e., to 
macroscopic polarization and strain) by a suitable change of the composition or by applying 
an external electric field; b) the existence of a frequency dependent slim hysteresis loop 
even above the transition temperature; c) the difference between field-cooled and zero field- 
cooled dielectric susceptibility, etc. All these findings indicate that the observed properties 
of relaxor ferroelectrics are nonequilibrium properties. Relaxor ferroelectrics thus represent 
a new low temperature state of polar dielectrics compared with conventional ferroelectrics. 

Experiments show that there is symmetry breaking on a nanometer scale as observed from 
Raman scattering ||^ and X-ray and neutron diffraction . This observation implies that 
polar clusters exist even well above and the observed properties of relaxors are strongly 
affected by the reorientations of the clusters. In particular, cluster reorientation in the 
applied electric field induces strong polar-strain coupling which makes relaxor ferroelectrics 



the candidate materials for the next generation of ultrasonic transducers ||16 



It has been proposed 0, in analogy with the cluster model of spin glasses, that the 
polar clusters behave like large superparaelectric dipole moments. The broad distribution 
of relaxation times for cluster orientations originates from the distribution of the potential 
barriers separating different orientational states. The superparaelectric model based on 
the assumption of independent clusters can not explain, however, the appearance of long 
range order with the change of the material composition or due to an applied external field. 
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The appearance of long range order also can not be explained by the inclusion of frustrated 
interactions between superparaelectric moments that lead to the dipole spin glass state. 
An alternative proposition PJT3|-p!5| relates the origin of relaxor behavior to the domain 
states induced by the static random field caused, for example, by charged composition fluc- 
tuations. Although this model seems very attractive, it encounters difficulties in explaining 
the observed cluster dynamics in the high temperature phase. 

We will show that the proposed model is capable to describe the anomalies of dielectric 
response of relaxor ferroelectrics and can be applied for systems possessing first or second 
order phase transitions or remaining only incipient ferroelectrics with very high dielectric 
constant and vanishing spontaneous polarization. 

The origin of polar clusters in relaxor ferroelectrics is still not understood. It has been 
proposed [|l^ that off-center ions might be responsible for relaxor properties in analogy 
with Ki-xLixTaO^ (KLT) where Li impurity ions occupy the off-center positions near the 
vacant K sites. It is well established jl^ that in the case of impurity induced relaxors like 



KLT the single impurity potential possesses multi-well structure which allows the thermal 
jumps of off-center ions between different potential minima. In contrast, in disordered 
complex perovskites the displacement of atoms from their centro-symmetrical positions is 
caused apparently by the charge compositional fluctuations, which violate the charge balance 
within the adjacent unit cells and lead to the additional electrostatic forces on the atoms. 
In this situation one would not expect that the single atom potential energy possesses multi- 
minima structure. Atom reorientations would more probably take place if they associate 
with the collective motion of atoms within small clusters, {i.e., the cluster potential energy 
is characterized by the multi-well (double well) structure). 

Such a picture is analogous to the situation in glasses where the origin of two level systems 
responsible for glass anomalies is associated with the double well potentials [0,0. It is 
generally assumed that the existence of the double well potentials is due to the disorder in 
glasses, so that local rearrangements of atoms might switch the system between two adjacent 
local energy minima. Attempts to detect double well potentials in computer simulations ||2^ 



based on the assumption of the motion of single atoms |2T|] have been unsuccessful. Instead 



the origin of double well potentials in glasses has been successfully explained ||2^ in terms 
of collective motion of atoms within small clusters. 

In this paper we develop further the cluster model of relaxor ferroelectrics taking ex- 
plicitly into account a) the broad distribution of local field experienced by each cluster due 
to cluster-cluster interaction and b) the broad distribution of potential barriers controlling 
cluster dynamics. We will show that the proposed model is capable of describing the anoma- 
lies of dielectric response of relaxor ferroelectrics and can be applied for systems possessing 
first or second order phase transitions or remaining an incipient ferroelectrics with very high 
dielectric constant and vanishing spontaneous polarization. 

II. MODEL 

Based on the analogy between disorder induced double well potentials in glasses and 
disordered relaxor ferroelectrics we will adopt the following picture of the cluster dynamics 
in relaxor ferroelectrics. Each minimum in a double well potential for a given cluster is 
characterized by the cluster dipole moment or cluster polarization Pd- The potential barrier 
between different minima determines the atom hopping frequency r^^. Clusters interact 
with each other, and interacting reorientable polarizable clusters should increase the crystal 
dielectric response. This conclusion is supported by the experiment in PST where the 
effect of disorder in the relative occupation of B sites by Sc or Ta atoms dramatically 
increases the dielectric constant. 

Time or frequency dependent dielectric response can be analyzed with the use of the 
master equation 

rIP 1 

^ = --(^d-P.?(i?L)), (1) 



well known from the theory of ferroelectrics of order- disorder type p3121], describing the 



relaxation of the polarization of each cluster to its quasi-equilibrium value P^{Ei) which 
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depends on the value of the local field El induced by other clusters at any moment of 
time. In general the local field is a time dependent random field. It also includes the 
contribution from the applied field E^x and the contribution from the static random fields 
caused by the material imperfections. 

In order to apply Eq.(|T]) to relaxor ferroelectrics one should perform the average over 
the distribution of relaxation times r and the distribution of the local fields El. For this 
purpose we rewrite Eq.(|lD in the equivalent integral form and take the average with respect 
to r, El and the initial cluster polarization Pd{0) to obtain 

Pit) = p(o)g(t) - dt'^^^^^p^^it'). (2) 

The function Q{t) = e"*/"^, where the overbar denotes the average over the relaxation time 
r, characterizes the slow nonexponential kinetics of the system. In addition 

P^'^{t) = j dEP:^{E)f{E,P{t)\ (3) 

where /(-E, P{t)) is the distribution function of the local field which depends parametrically 
on the value of the average polarization of the system P{t). We will consider below f{E, P) 
in the form f{E, P) = f{E — 7P — 7eo-Eex/4'/r). This form of / is consistent with the mean 
field approximation f{E, P) = 6{E — ■yP — ■yeoEex/4:7!'), where 6 is the delta-function and 7 is 
the local field phenomenological parameter. In disordered systems the effect of composition 
fluctuations leads to a deviation from the simple mean fleld picture that can be taken into 
account by the replacement of the (5-function by the function / with flnite width. The shape 
and the width of f{E) depend on the explicit form of the cluster-cluster interactions as well 
as on static random fields caused by the material imperfections. 

The value 7eo-Eea/47r is the local field induced by the external field in the dielectric 
media with the dielectric constant eo >> 1 |^ {i.e., we assume that the polar clusters 
are distributed in a highly polarizable dielectric media). In relaxor ferroelectrics which are 
mainly perovskite-based highly polarizable materials, the typical values are cq > 10^. 

Eq-dl) can be applied for the analysis of different experimental situations {e.g., decay 
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of the polarization, the difference between field-cooled and zero field-cooled dielectric sus- 
ceptibility, the effect of a frequency dependent hysteresis loop, etc). In this paper we will 
concentrate, for illustration, on the calculation of the frequency dependent linear dielectric 
susceptibility. In order to calculate the steady state susceptibility in the presence of an ad- 
ditional time independent field Ej°J we write E^^{t) = E'fJ + ^^^e^"^* and P{t) = Ps + Pi(t), 
where Pg is the time independent polarization induced by the field E^^ (for E^^ = 0, we de- 
note Ps as the spontaneous polarization). We obtain from Eq.(^ the following self consistent 
equation for Pg 

Ps = jdEP:,\E)f{E,Pg). (4) 

Assuming P(0) = P, in Eq.(|D, considering a linear expansion of P^'?(t) with respect to 
and Pi{t), and taking Laplace transform of Eq.(0), we obtain 

1 - k{T)Q{uj,T) 

where we used the definition of the dielectric constant e{uj,T) = ATxdPi{uj) / dE^J + eo- In 
Eq.® 

K(T) = J dEP:m^J^§^, (6) 

and Q{uj,T) is the Laplace transform of dQ{t)/dt given by 



Q(u;,T) = (1/(1 +z^r)) (7) 
Assuming an Arrhenius or Vogel-Fulcher (VF) law for r 

r(T) = roexp[f//(T-To)], (8) 
we may calculate Q{u,T). Thus, we obtain for Q'{lj,T) p6| 

/•{T-To)ln^ 

Q\u;,T)= ^^'dUgiU), (9) 







where g(U) is the distribution function of the potential barriers. 



In the spirit of Landau phenomenological theory one can further expand the right hand 
side of Eqs. (^) and in a power series with respect to assuming that the apphed field 
is sufficiently small, and present k as 

K = ai + 3asP^ + ba^P^ (10) 

Eqs-dl), (H), and (|T0]) give the phenomenological description of the dynamical response of 
relaxor ferroelectrics and can be applied for systems possessing first or second order phase 
transitions or remaining as only incipient ferroelectrics with very high dielectric constant 
and vanishing spontaneous polarization. 

III. AMOMALIES IN THE DIELECTRIC RESPONSE 

In this section we illustrate a capability of the model in the description of the dielectric 
response of relaxor ferroelectrics. 

Relaxors with incipient ferroelectric order. In order to simultaneously reproduce the high 
values of the dielectric constant in relaxor ferroelectrics and the absence of spontaneous 
polarization (like that in PMN or PST with vacancies) one should assume that 03 < and 
ai (T) — > 1 remaining, however, less than 1 at all temperatures. We chose for the illustrative 
calculations ai(T) = 0.95 tanh(0.5/T) (which, according to Eq.(||), reproduces Curie-like 
high temperature behavior of e(0, T) and its high saturation value at low temperatures) and 
g{U) = (0.5 + f/^)^. We assume also the Arrhenius temperature dependence of t(T), 

i.e. To = in Eq.(|D. Using the chosen values of ai(T), g{U), and Tq we calculated the 
dielectric permittivity e as a function of temperature for different values of In(^). 

The calculated real part of permittivity is shown in Fig.l Note that the frequency depen- 
dent maximum of e' is not just a relaxation maximum. It originates from the competitive 
temperature dependences of Q'{T) and ai(T). The behavior obtained for e' is in reasonable 
qualitative agreement with the experiment For a more detailed comparison with 

the experiment one needs to find the functions Q{uj^ T) and cii(T) corresponding to the best 
fit with the experimental data. 
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Manifestation of the first order phase transition. In some relaxor ferroelectrics (e.g., 
disordered PST or KLT above the critical concentration of Li ions) there is evidence of a 
first order phase transition, which manifests itself in the sharp drop of the dielectric constant 
at the temperature below the position of the relaxation maximum. In order to reproduce 
this behavior within the proposed phenomenological theory one should assume 03 > 0. For 
example for 03 = 1.7, 05 = —11 and the values of ai(T) and Q' being the same as used above, 
the first order phase transition occurs at Tc ~ 0.2 (in the chosen dimensionless units). 

The calculated values of e' are shown in Fig. 2. A remarkable feature of the e' temperature 
dependence is that the relaxation maximum of e' approaches the phase transition temper- 
ature with a decrease of the frequency until it finally disappears transforming to a sharp 
peak like that in conventional ferroelectrics. This behavior is in qualitative agreement with 
that observed in PST. 

Reconstruction of the relaxation function Q{uj,T). PST with B-site chemical disorder 
undergoes a first order relaxor-ferroelectric phase transition at Tc ~ 269K ^ . The dielectric 
permittivity of PST shows pronounced frequency dispersion with the position of a frequency 
dependent maximum obeying the VF law |]^ c<j = u!oexp[—U/(Tm — Tq)], where u is the 
frequency of the applied field and is the temperature of the permittivity maximum 
corresponding to the frequency u . It has been widely accepted (see, e.g., |^^) that the VF 



type relation for the permittivity maximum is a consequence of the VF law (P) for t(T). This 
assumption has been argued recently by Tagantsev ||2^ who proposed that the observation in 



PST of the VF frequency dependence of can be explained with the use of the Arrhenius 
temperature dependence of r taking into account the existence of the first order phase 
transition and the fact that in PST Tq ^ Tc. An important conclusion of Tagantsev's 
analysis is the indication that might be infiuenced significantly by the values of static 
permittivity, not only by the specific temperature dependence of r. However, in order to be 
more conclusive one should analyze the shape of the frequency dependent permittivity as a 
function of temperature for different frequencies, not just the frequency dependence of the 
permittivity maximum. Such an analysis is given below with the use of Eq.(^). 
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First, we extrapolated the experimental data |^ and extracted the static permittivity 
above as shown in Fig. 3 (dotted curve). Then, using Eq.(||) and the fact that at u; = Q' 
is equal to 1, we obtained the values of /«(T) above Tc. We used the value eo ~ 500 for the 
host lattice permittivity based on the experimental low temperature values of e(T) where 
the polar clusters are frozen. With the values obtained for k(T) and the experimental values 
of the frequency dependent permittivity e'{uJi,T) at frequency ui = lOKHz (curve 1) we 
have calculated the values of Q'{u!i,T) using Eq.(||). According to Eq.(|D Q'{u,T) is a 
function of (T — Tq) In(co'ro) that results in the scaling relation Q'{u!,T) = Q'{ui,Ti) with 
Ti = (T — To)^^l^^ + Tq. Using this relation one can reconstruct the values of Q'{uj, T) at 
other frequencies and, therefore, the values of e'(co',T). The results of such a reconstruction 
are presented in Fig.3 for the upper lOOHz (curve 2) and lower IMHz (curve 3) boundary 
frequencies used in the experiment. We used tq ~ 10~^^ in accordance with Ref. [^. 

One can see from Fig.3 that the approach reproduces rather well the shape of the tem- 
perature dependence of the permittivity at different frequencies of the applied field. The 
parameter of the fit Tq is found to be Tq ~ 2b%K which is in reasonable agreement with the 
value To ~ 2QbK estimated in Ref. ||^. This result means that the employment of the VF 
law for t{T) is very crucial for obtaining the correct temperature and frequency dependence 
of the permittivity for PST. The VF relation for t(T) results simultaneously in the very 
fast temperature dependence of the permittivity and its rather slow dependence on ln(ti;ro), 
observed in the experiment. At the same time an Arrhenius like dependence of t(T) would 
result in an extremely rapid dependence of the permittivity on In(co'ro), which is inconsistent 
with the experimental data. 

We are grateful to A. Khachaturian, A. Tagantsev, and J. Toulouse for the useful discus- 
sions. This work is supported by the National Science Foundation and the Army Research 
Office. 
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FIGURES 



FIG. 1. Real part of the dielectric susceptibility in relaxors with incipient ferroelectric order. 

FIG. 2. Effect of the first order phase transition on the dielectric response in relaxor ferro- 
electrics. 

FIG. 3. e'(T) for disordered PST above Tc for different frequencies. The vertical line shows 
the phase transition temperature. Solid lines 2,3 and the circles are respectively the reconstructed 
and the experimental values of e'{T). The error bars indicate the uncertainties of the reconstruc- 
tion due to the uncertainties in the values obtained for k{T) by the extrapolation(dotted line) of 
experimental data. 
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